A Radio Galaxy at z = 5.19 1 



Wil van Breugel 2 , Carlos De Breuck 2 ' 3 , S.A. Stanford 2 , 
Daniel Stern 4 , Huub Rottgering 3 , & George Miley 3 



ABSTRACT 

We report the discovery of the most distant known AGN, the radio galaxy 
TN J0924— 2201 at z = 5.19. The radio source was selected from a new sample 
of ultra-steep spectrum (USS) sources, has an extreme radio spectral index 
a 365MHz = —1-63, and is identified at near-IR wavelengths with a very faint, 
K = 21.3 ± 0.3 object. Spectroscopic observations show a single emission line 
at A ~ 7530 A which we identify as Lya. The AT-band image, sampling rest 
frame [/-band, shows a multi-component, radio-aligned morphology, typical 
of lower-redshift radio galaxies. TN J0924— 2201 extends the near-IR Hubble, 
or K — z, relation for powerful radio galaxies to z > 5, and is consistent with 
models of massive galaxies forming at even higher redshifts. 



Subject headings: cosmology: early universe - galaxies: active - galaxies: 
distances and redshifts - galaxies: formation - galaxies: individual: TN 
J0924-2201 



1. Introduction 

How did the first objects form after the Big Bang? In hierarchical cosmogonies 
(e.g., Turner 1998), the first gravitationally bound systems may have been stars and small 
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star-forming systems which merge to form galaxies in large dark matter halos. Arising 
from the end products of stellar evolution and mergers, central black holes could grow to 
become extremely massive. However, it is not clear how this process would work at very 
high redshifts, where little time is available. It has been suggested that primordial black 
holes may form well before their host galaxies (Loeb 1993). In any case, accretion events 
fueling massive black holes are thought to manifest themselves as active galactic nuclei 
(AGN; e.g., Rees 1984). Due to their extreme luminosity, AGN are convenient beacons for 
exploring these formative, 'Dark Ages' of our Universe. 

Extragalactic radio sources have played an important role in identifying active galaxies 
at high redshifts. The most distant known galaxies have consistently been radio-selected 
until only very recently. In this Letter we report the discovery of a radio galaxy at z = 5.19. 
At this redshift it is the most distant known AGN, surpassing even quasars for the first 
time in 36 years. Throughout this paper we use H = OS/i^kuts^ 1 Mpc -1 , f2 M = 0.3, and 
A = 0. For these parameters, 1" subtends 7.0 h$ 5 kpc at z — 5.19 and the Universe is only 
1.08 Gyr old, corresponding to a lookback time 91.1% of the age of the Universe. 



2. Source Selection 

The most efficient method to find high-redshift radio galaxies (HzRGs) is to combine 
two well-known techniques. The first is to select radio sources with ultra-steep spectra 
(USS) at radio wavelengths, i.e. very red radio colors (e.g., Chambers, Miley, & van Breugel 
1990). Most powerful radio galaxies have radio spectral energy distributions which steepen 
with frequency. Therefore, at fixed observing frequencies more distant sources exhibit 
steeper spectra (e.g., van Breugel et al. 1999). 

A second selection criterion relies upon the magnitude-redshift relationship at infrared 
wavelengths, or K — z Hubble diagram, for powerful radio galaxies (Figure [I]). At low 
redshifts (z < 1), powerful radio galaxies are uniquely associated with massive galaxies. 
The well-behaved K — z diagram suggests that such galaxies can be found through near-IR 
identification. This has been confirmed by the discovery of many 3 < z < 4.4 radio galaxies 
which approximately follow the K — z relationship, even to the highest redshifts and despite 
significant morphological evolution (van Breugel et al. 1998). 

Using several new, large radio surveys we constructed a USS sample 
(S„ cx z/^o4S z < -1.30; De Breuck et al. 1999 [DB99]) which is much larger, 
more accurate, and reaches fainter flux density limits than previous such samples. 
TN J0924— 2201, with a^j^ = —1.63 ± 0.08, is among the steepest sources of our sample. 
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VLA observations at 4.85 GHz show the source is a slightly resolved 1'.'2 double, with 
Sa.85GHz = 8.6 ±0.5 mJy, centered at a J20 oo = 09 h 24 m 19 f 592, 5 J20 oo = -22°01'41'. / 5 (Figure |). 

3. Observations 

We obtained K s images of TN J0924-2201 using NIRC (Matthews & Soifer 1994) at 
the Keck I telescope. We integrated for 32 minutes on UT 1998 April 18 in photometric 
conditions with 0'.'5 seeing, and again for 32 minutes on UT 1998 April 19 through light 
cirrus with 0'.'6 seeing. The observing procedures, calibration, and data reduction techniques 
were similar to those described in van Breugel et al. (1998). The final image comprising 
3840 s of on-source integration is shown in Figure [| Using circular apertures of 2 / .'l 
diameter, encompassing the entire object, we measure K = 21.15 for night 1, and 21.45 for 
night 2. We estimate that K = 21.3 ± 0.3. If TN J0924-2201 is at z = 5.19 (§4), then 
redshifted [O II] A 3727 at A = 2.307/mi would be included in the K s passband and some of 
the X-band flux might be due to line emission. 

We obtained spectra of TN J0924-2201 through a 1"5 wide, 3' long slit using LRIS 
(Oke et al. 1995) at the Keck II telescope. The integration times were 5400 s on UT 1998 
December 19 (position angle 0°) and 4400 s on UT 1998 December 20 (position angle 180°); 
both nights were photometric with / .'6 seeing. The observations used the 150 lines mm -1 
grating (Abiazc ~ 7500 A; AAfwhm ~ 17 A), sampling the wavelength range 4000 A to 
1/im. Between each 1800 s exposure, we reacquired offset star A (see Fig. 2), performed 
20" spatial shifts to facilitate removal of fringing in the reddest regions of the spectra, and 
blind offset the telescope to return TN J0924— 2201 within the slit. We calculated the 
dispersion using a NeAr lamp spectrum taken immediately subsequent to the observations 
(RMS variations of 0.50 A), and adjusted the zero point according to telluric emission 
lines. Final wavelength calibration is accurate to 1 A. The spectra were flux calibrated 
using observations of Feige 67 and Feige 110 obtained on each night and were corrected for 
foreground Galactic extinction using a reddening of Eb-v — 0.0168 determined from the 
dust maps of Schlegel, Finkbeiner, & Davis (1998). 

We find a strong, single emission line at A ~ 7530 A which shifts by ~ 16 A between 
the two nights. (Figure |3|; Table 1). The cause of the line offset is unclear, though it 
may be related to problems LRIS was experiencing with slippage in the movable guider at 
the time of the observations. The relative brightnesses of other sources on the slit vary 
between each 1800 s observation, indicating that despite our precautions of reacquiring the 
target after each exposure, guider slippage must have caused some variations in telescope 
offsetting. These slight pointing changes may have caused the slit to sample different 
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regions of spatially-extended, line-emitting gas. Indeed, TN J0924— 2201 shows two 
separate components at K (Figure |2[), and emission-line regions of HzRGs are known to be 
kinematically complex (Chambers, Miley & van Breugel 1990; van Ojik et al. 1997). 

Line parameters are measured with a Gaussian fit to the emission line and a flat (in 

Fx) fit to the continuum (Table 1). Equivalent width values were derived from a Monte 

Carlo analysis using the measured line flux and continuum values with errors, subject 

to the constraint that both are positive. For UT 1998 Dec. 19, when no continuum was 

reliably detected, we quote the 90% confidence limit, W ohs > 2760 A. For UT 1998 Dec. 

20, when continuum was marginally detected, we quote the 90% confidence interval, 
W obs = 71Q _ 155Q ^ 



4. Redshift Determination 

As discussed by Dey et al. (1998) and Weymann et al. (1998) for two z > 5 Lya-emitting 
field galaxies, a solitary, faint emission line at red wavelengths is most likely to be either 
low-redshift [O II] A 3727 or high-redshift Lya. Similar arguments are even more persuasive 
for HzRGs because of their strong, rich emission line spectra. For example, if the line at ~ 
7530 A were [O II] at z — 1.020 then composite radio galaxy spectra (McCarthy 1993; Stern 
et al. 1999a) indicate that the TN J0924-2201 spectrum should have shown C II] A 2326 
at 4699 A with w 40 - 70% the strength of [O II], and Mg II A 2800 at 5653 A with 
~ 20 — 60% the strength of [O II]. Similar arguments rule out identifying the emission line 
with Ha at z — 0.147 or [O III] A 5007 at z = 0.504, since in these cases even stronger 
confirming lines should have been seen. 

The large equivalent widths also argue against identifying the emission line with [O II] 
at z = 1.020, implying Wffi > 1370 A (night 1) and 350 < Wffi < 770 A (night 2). Radio 
galaxy composites typically have rest-frame [O II] equivalent widths of ~ 130 A (McCarthy 
1993; Stern et al. 1999a), though active galaxies with extreme W^J] are occasionally 
observed (W^% ~ 750 A ; Stern et aL 1999b). The equivalent width of TN J0924-2201 is 
more typical of high-redshift Lya which is often observed with rest frame values of several 
x 100 A in HzRGs (Table 2). We also note that the observations from the second night 
show that Lya is attenuated on the blue side, presumably due to associated and intervening 
hydrogen gas, as is commonly observed in HzRGs (e.g., van Ojik et al. 1997; Dey 1997) 
and normal star- forming galaxies at z > 5 (e.g., Dey et al. 1998). 

Finally, the faint K-band magnitude of TN J0924— 2201 conforms to the extrapolation 
of the K — z relation to z > 5 (Figure [I]). Identifying the emission line with [O II] would 
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imply a severely underluminous HzRG (by 3-4 mag). Therefore, the most plausible 
identification of the emission line in TN J0924— 2201 is with Lya at a (mean) observed 
wavelength of 7530 A and z = 5.19. Table 1 gives the dereddened emission-line fluxes. 

5. Discussion 

Among all known z <; 3.8 HzRGs, TN J0924— 2201 is fairly typical in radio luminosity, 
equivalent width, and velocity width (Table 2). But this source has the steepest radio 
spectrum, consistent with the a — z relationship for radio galaxies (e.g., Rottgering et 
al. 1997). TN J0924— 2201 also has the smallest linear size, perhaps indicating that the 
source is relatively young and/or embedded in a denser environment compared to the other 
HzRGs, commensurate with its large velocity width (van Ojiket al. 1997) and very high 
redshift. Together with 8C 1435+63, TN J0924— 2201 appears underluminous in Lya, which 
might be caused by absorption in a relatively dense cold and dusty medium. Evidence 
for cold gas and dust in some of the most distant HzRGs has been found from sub-mm 
continuum and CO-line observations of 8C 1435+63 and 4C 41.17 (e.g., Ivison et al. 1998). 

Our observations of TN J0924— 2201 extend the Hubble K — z diagram for powerful 
radio galaxies to z = 5.19. Simple stellar evolution models are shown in Figure [l] for 
comparison with the HzRG. Despite the enormous /c-correction effect (from U Tes t at z — 5.19 
to -ft^est at z = 0) and strong morphological evolution (from radio-aligned to elliptical 
structures) , the K — z diagram remains a powerful phenomeno logical tool for finding radio 
galaxies at extremely high redshifts. Deviations from the K — z relationship may exist 
(Eales et al. 1997; but see McCarthy 1998), and scatter in the K — z values appears to 
increase with redshift. 

The clumpy, radio-aligned U rest morphology resembles that of other HzRGs (van 
Breugel et al. 1998; Pentericci et al. 1998). If the continuum is dominated by star light, as 
appears to be the case in the radio-aligned HzRG 4C 41.17 at z — 3.798 (Dey et al. 1997), 
then M(U) = -24.4 for TN J0924-220lthen we can derive a SFR of ~200 M yr" 1 , 
assuming a Bruzual & Chariot (1999) GISSEL stellar evolution model with metallicity 
Z = 0.008, no extinction, and a Salpeter IMF. This SFR value is highly uncertain due to 
the unknown, but competing, effects of extinction and [O II] emission-line contamination, 
but is not unreasonable. It is 2.5 times less than in 4C 41.17, which has M(U) = —25.2 
using the same aperture (Chambers et al. 1990). TN J0924— 2201 may be a massive, active 
galaxy in its formative stage, in which the SFR is boosted by induced star formation 
(e.g., Dey et al. 1997). For comparison other, 'normal' star forming galaxies at z > 5 have 
10 - 30 times lower SFR (~ 6 - 2OM ?/r~ 1 ; Dey et al. 1998; Weymann et al. 1998; Spinrad 
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et al. 1998). 

Recent z ~ 3 and z ~ 4 Lyman-break galaxy observations have suggested a possible 
divergence of star formation and AGN activity at high redshift (Steidel et al. 1999), contrary 
to what was previously thought (e.g., Haehnelt, Natarajan & Rees 1998). However, if 
starbursts and AGN are closely coupled, as suggested to explain the ultraluminous infrared 
galaxies (Sanders & Mirabel 1996), then young AGN may inhabit especially dusty, obscured 
galaxy systems. To obtain a proper census of the AGN population at the very highest 
redshifts therefore requires samples which avoid optical photometric selection and extinction 
bias, such as our cm-wavelength/A'-band radio galaxy sample. 

As emphasized by Loeb (1993), if massive black holes form in a hierarchical fashion 
together with their host galaxies, this process must be quick and efficient, as available 
timescales are short: at z = 5.19 the Universe is only 1 Gyr old. It is unclear how this could 
be done, so other models, where primordial massive black holes form soon after the Big 
Bang and prior to the beginning of galaxy formation, may require additional investigation. 

We thank G. Puniwai, W. Wack, R. Goodrich and R. Campbell for their expert 
assistance during our observing runs at the W.M. Keck Observatory, and A. Dey, J.R. 
Graham and H. Spinrad for useful discussions. The work by W.v.B., C.D.B. and S.A.S. 
at IGPP/LLNL was performed under the auspices of the US Department of Energy under 
contract W-7405-ENG-48. W.v.B. also acknowledges support from NASA grant GO 5940, 
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Table I. Spectroscopic Measurements of TN J0924-2201 



UT Date 1998 Dec 19 1998 Dec 20 





5.202 ± 0.002 


5.188 ±0.001 


A [A] 


7539.2 ± 1.9 


7522.9 ± 1.5 


F Lya [10- 17 erg cm- 2 s" 1 ] 


3.4 ±0.6 


3.5 ±0.6 


F cont [ 10 -2l crg cm -2 g -l 


3.3 ±5.7 


33.0 ±4.9 


w°^ s a [A] 


> 2760 


710 - 1550 


FWHM Lya [km s- 1 ] 


1574 ± 192 


1474 ± 106 



Table 2. Physical parameters of HzRGs 



Name 


z 


L Lya x 10~ 43 
(crg s -1 ) 


isMMEz X 10 36 
(org s -1 ) 


1.4GHz 
a 365MHz 


Lye* 

(A) 


FWHM LyQ 
(km s- 1 ) 


Radio size 
(kpc) 


Reference 


TN J0924-2201 


5.19 


1.3 


7.5 


-1.63 


>115 


1500 


8 


This paper; DB99 


6C 0140+326 


4.41 


11 


1.3 


-1.15 




1500 


19 


Rawlings et al. (1996) 


8C 1435+63 


4.25 


3.2 


11 


-1.31 


670: 


1800 


28 


Spinrad et al. (1995) 


TN J1338-1942 


4.13 


20 


2.3 


-1.31 


>700: 


1300 


37 


DB99 


4C 41.17 


3.798 


12 


3.3 


-1.25 


100 


1400 


99 


Dey et al. (1997) 


4C 60.07 


3.79 


16 


4.1 


-1.48 


150 


2900: 


65 


Rottgering et al. (1997) 



12 



10 
0.01 




3C sample 
6C sample 

NIRC HzRG 
1 Gyr burst @ z f =20 
0.1 Gyr burst @ z f =20 



0.1 



10 



Redshift z 



Fig. 1. — Hubble K — z diagram. Filled triangles are Keck measurements of HzRGs from 
van Breugel et al. (1998), the large triangle is TN J0924— 2201, and all other photometry 
is from Eales et al. (1997). Magnitudes are aperture-corrected to a 64 kpc metric diameter 
using H = GSkms- 1 Mpc" 1 and = 0.30, for which TN J0924-2201 has K = 20.7 ± 0.3. 
Two stellar evolution models from Bruzual & Chariot (1999), normalized at z < 0.1, are 
plotted, assuming parameters as shown. 
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Fig. 2. — Keck/NIRC AT-band image of TN J0924— 2201, with radio contours superposed. 
Star 'A', offset from the HzRG by Aa = 7"2 (E) and A5 = 1877 (N), was used for blind- 
offsettine for the SDectrosconic observations. 
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Fig. 3.— Keck spectra of TN J0924-2201. Top panel is night 1, bottom is night 2. Spectra 
were extracted using a 1'.'5 x 1"5 aperture. Inserts illustrate the emission line morphology, 
which is attenuated on the blue edge for night 2. The predicted wavelength of N V A 1240 
emission is indicated. 



